Measurement of the masses and widths of the bottom baryons and 

T. Aaltonen,^^ B. Alvarez Gonzalez^ ,^ S. Amerio,^" D. Amidei,^^ A. Anastassov^/^ A. Annovi,^'' J. Antos,^^ 
G. Apollinari/^ J. A. Appel,^^ T. Arisawa,^"* A. Artikov/^ J. Asaadi,''^ W. Ashmanskas,^^ B. Auerbach,^"^ 
A. Aurisano.'^s F. Azfar,^^ W. Badgett/^ T. Bae,^^ A. Barbaro-Galtieri,^^ V.E. Barnes,''^ B.A. Barnett,^^ 
P. Barria'''',42 p. Bartos/^ M. Bauce-'^-'' F. Bedeschi,^^ Behari,23 G. Bellettinifs ,^^2 Bellinger, 
D. Benjamin/-* A. Beretvas/^ A. Bhatti,*^ D. Bisello^-'',''^ I. Bizjak,^^ K.R. Bland,^ B. Blumenfeld,^^ A. Bocci," 

A. Bodek,'*^ D. Bortoletto,-*-* J. Boudreau,'*^ A. Boveia," L. Brigliadori^^^ (j Bromberg,^^ E. Brucken,^! 
J. Budagov,i3 H.S. Budd,*"^ K. Burkett,*^ G. Busetto^-'',*'^ p. Bussey,!^ A. Buzatu,^! A. Calamba,i° C. Calancha,^ 
S. Camarda," M. Campanelli,^^ M. Campbell,^^ Canelli",*^ B. Carls,22 D. Carlsmith.^*^ R. Carosi,*^ 
S. CarriUo™,!*^ S. Carron,!^ B. Casal^^^^ M. Casarsa,50 A. Castro'=%6 P. Catastini,^" D. Cauz,^^ V. Cavaliere,22 
M. Cavalli-Sforza,4 A. Cerri^,26 L. Cerrito^28 Y.C. Chen,i M. ChertokJ G. Chiarelli,*^ q Chlachidze,!^ 
F. Chlebana/^ K. Clio,^^ D. Chokheli," W.H. Cliung,^^ y.S. Chung,^^ M.A. Ciocci'''',''^ A. Clark,!^ C. Clark 
G. Compostella-''^,40 M.E. Convery/^ J. Conway/ M.Corbo/^ M. Cordelli/^ C.A. Cox/ D.J. Cox/ F. Crescioli 
J. Cuevas^9 R. Culbertson/'^ D. Dagenhart/'^ N. d'Ascenzo™/^ M. Datta/^ P. de Barbaro/'^ M. DeirOrso^s 
L. Demortier/6 M. Deninno/ F. Devoto/^ M. d'Errico^^/" A. Di Canto^f /^ b. Di Ruzza/^ J.R. Dittmann 
M. n'Onofrio 27 Q. nona,ti99.42 P. Donsr M. noritro.S" T. Doriffo.^O K. Ebina.S* A. Elamn>9 A. Enniff 32 



J. Uuevas"-, K. Uuibertson,"" U. UagenJiart/" IN. a Ascenzo-"/" M. Uatta, " ae tiarbaro, " M. Uell Urso''", ' 
L. Demortier/6 M. Deninno/ F. Devoto/^ M. d'Errico^^/" A. Di Canto^f /2 b. Di Ruzza/^ J.R. Dittmann/ 
M. D'Onofrio/^ S. Donati^f /2 Dong/^ M. Dorigo/" T. Dorigo/° K. Ebina/"^ A. Elagin/^ A. Eppig/2 
R. Erbacher/ S. Errede/2 N. Ershaidat'^''/^ R. Eusebi/^ S. Farrington/^ M. Feindt/'' J.P. Fernandez/^ 
R. Field/6 G. Flanagan"/^ R. Forrest/ M.J. Frank/ M. Franklin/" J.C. Freeman/^ Y. Funakoshi/* I. Furic/*^ 

M. Gallinaro/6 j.E. Garcia/* A.F. Garfinkel/* P. Garosi''''/2 jj. Gerberich/2 E. Gerchtein/'^ S. Giagu/^ 
V. Giakoumopoulou/ P. Giannetti/2 K. Gibson/^ C.M. Ginsburg/^ N. Giokaris/ P. Giromini/^ G. Giurgiu/^ 

ir 13 t-v m -.^.i-ISa^ /-i„ij35t-. 49at /-i„ij„„i :jil6 A /-i„i 15 ri 9 



V. Giakoumopoulou/ P. Giannetti/2 Gibson/^ C.M. Ginsburg/^ N. Giokaris/ P. Giromini/^ G. Giurgiu/^ 
V. Glagolev/3 D. Glenzinski/^ M. Gold/^ D. Goldin/^ N. Goldschmidt/^ A. Golossanov/^ G. Gomez/ 
G. Gomez-Ceballos/° M. Goncharov/o O. Gonzalez/^ I. Gorelov/^ A.T. Goshaw/* K. Goulianos/^ S. Grinstein^ 
C. Grosso-Pilcher/i R.C. Group^^ J. Guimaraes da Costa/o S.R. Hahn/^ E. Halkiadakis/® A. Hamaguchi/^ 
J.Y. Han/5 F. Happacher/^ K. Hara/* D. Hare/** M. Hare/2 R.F. Harr/^ K. Hatakeyama/ C. Hays/^ 
M. Heck/4 J. Heinrich/i M. Herndon/^ S. Hewamanage/ A. Hocker/^ W. Hopkins'?/^ D. Horn/-* S. Hon/ 
R.E. Hughes/'^ M. Hurwitz/^ U. Husemann/'^ N. Hussain/* M. Hussein/^ J. Huston/^ G. Introzzi/2 M. lori" 
A. IvanovP/ E. James/^ D. Jang/° B. Jayatilaka/* E.J. Jeon/^ S. Jindariani/^ M. Jones/'' K.K. Joo/^ 
S.Y. Jun/o T.R. Junk/^ T. Kamon25/9 p p Karchin/^ A. Kasmi/ Y. Kato"/^ W. Ketchum/i J. Keung/* 
V. Khotilovich/9 B. Kilminster/^ D.H. Kim/^ H.S. Kim/^ J.E. Kim/^ M.J. Kim/^ S.B. Kim/^ S.H. Kim/* 
Y.K. Kim/i Y.J. Kim/5 ^ Kimura/* M. Kirby/^ S. Klimenko/^ K. Knoepfel/^ K. Kondo*/^ D.J. Kong/^ 
J. Konigsberg/6 A.V. Kotwal/'' M. Kreps/* J. Kroll/^ D. Krop/* M. Kruse/^ V. Krutelyov/*^ T. Kuhr/* 

Kuratfl 51 Q Kwa.ncr ^-'^ A T T.aafia.npn T, Labaro-a"^ 29 Q T,a.mi '*2 Q T.ammpl TVT T.anraster 28 VI T, T.anHp 



J. Konigsberg/6 A.V. Kotwal/'' M. Kreps/'' J. Kroll/^ D. Krop/* M. Kruse/^ V. Krutelyov/*^ T. Kuhr/^ 
Kurata/* S. Kwang/-' A.T. Laasanen/* L. Labarga"^/^ S. Lami/2 S. Lammel/^ M. Lancaster/^ R.L. Lander/ 
Lannon^'/e A. Lath/^ G. Latino''''/2 T. LeCompte/ E. Lee/^ H.S. Lee?/* J.S. Lee/^ S.W. Lee''^*9 S. Leo9f/2 
S. Leone/2 J.D. Lewis/^ A. Limosani*/-* C.-J. Lin/^ M. Lindgren/^ E. Lipeles/* A. Lister/^ D.O. Litvintsev/^ 
C. Liu/3 H. Liu/3 44 p Lj^^i5 Lockwitz/^ A. Loginov/^ D. Lucchesi^^/° J. Lueck/* P. Lujan/^ 

P. Lukens/5 G. Lungu/^ J. Lys/^ R. Lysak/i2 Madrak/^ K. Maeshima/^ P. Maestro'*'' /2 S. Malik/^ 
G. Manca"/'' A. Manousakis-Katsikakis/ F. Margaroli/'^ C. Marino/"* M. Martinez/ P. Mastrandrea/'^ 
K. Matera/2 M.E. Mattson/^ A. Mazzacane/^ P. Mazzanti/ K.S. McFarland/^ P. Mclntyre/^ R. McNulty^/^ 
A. Mehta/'^ P. Mehtala/* C. Mesropian/*^ T. Miao/^ D. Mietlicki/2 A. Mitra/ H. Miyake/* S. Moed/^ 
N. Moggi/ M.N. Mondragon^/s C.S. Moon/5 R. Moore/5 M.J. Morello"/2 J. Morlock/"* P. Movilla Fernandez/5 
A. Mukherjee/5 Th. MuUer/* P. Murat/5 M. Mussini'=^6 J. Nachtman"/5 Y. Nagai/* J. Naganoma/"* L Nakano/'^ 
A. Napier/2 J. Nett/^ C. Neu/3 M.S. Neubauer/2 J. Nielsen'' /^ L. Nodulman/ S.Y. Noh/5 0. Norniella/2 
L. Oakes/9 S.H. Oh/* Y.D. Oh/5 j oksuzian/3 T. Okusawa/^ R. Orava/* L. Ortolan/ S. Pagan Griso^-'^/o 
C. Pagliarone/° E. Palencia^/ V. Papadimitriou/5 A. A. Paramonov/ J. Patrick/5 q Pauletta'''' /° M. Paulini/° 
C. Paus,3° D.E. Pellett/ A. Penzo/" T.J. Phillips/* G. Piacentino/2 E. Pianori/* J. Pilot/^ K. Pitts/2 C. Plager/ 
L. Pondrom/6 S. Poprockif/5 K. Potamianos/* F. Prokoshin=/*3 Pranko/^ F. Ptohos''/^ G. Punzi9s/2 
A. Rahaman/3 V. Ramakrishnan/^ N. Ranjan/* L Redondo/^ P. Renton/^ M. Rescigno/'' T. Riddick/^ 
F. Rimondi'=^^ L. Ristori*2/5 A. Robson/^ T. Rodrigo/ T. Rodriguez/* E. Rogers/2 S. RoUi/52 Roser/5 
F. Ruffini''''/2 A. Ruiz/ J. Russ/° V. Rusu/5 A. Safonov/^ W.K. Sakumoto/5 Y. Sakurai/* L. Santi'='=/° 
K. Sato/* V. Savehev'^/5 A. Savoy-Navarro°°/5 P. Schlabach/5 A. Schmidt/* E.E. Schmidt/5 T. Schwarz/5 
L. Scodellaro/ A. Scribano''''/2 F. Scuri/2 S. Seidel/5 Y. Seiya/** A. Semenov/3 F. Sforza''''/2 S.Z. Shalhout/ 
T. Shears/'^ P.F. Shepard/3 M. Shimojima"/* M. Shochet/* L Shreyber-Tecker/* A. Simonenko/3 P. Sinervo/* 



2 



K. Sliwa,52 J Smith/ F.D. Snider,!^ A. Soha,!^ V. Sorin,'* H. Song,43 p. Squillacioti'''\42 M. Stancari/^ 
R. St. Denis,!^ B. Stelzer,^! O. Stelzer-Chilton,^! D. Stentz^/^ j strologas,^^ G.L. Strycker,32 Y. Sudo,^! 
A. Sukhanov/^ I. Suslov/^ K. Takemasa,^^ Y. Takeuchi,^! J. Tang," M. Tecchio,^^ P.K. Teng/ J. Thoms/^ 
J. Thome,io G.A. Thompson,22 E. Thomson,^! D. Toback,49 S. Tokar,!^ K. Tollcfson,-'^^' T. Tomura,^! D. ToncUi,!"^ 
S. Torre/'^ D. Torretta,!^ P. Totaro,'^° M. Trovato",'^^ p_ Ukegawa,^! S. Uozumi,^^ A. Varganov,^^ F. Vazquez", 
G. Velev,i5 C. Vellidis,^^ ^ vidal,*^ I. Vila,^ R. Vilar,^ J. Vizan,^ M. Vogel,^^ q Volpi,i^ R Wagner,^ 
R.L. Wagner,i5 T. Wakisaka,^^ R. Wallny,^ S.M. Wang,i A. Warburton.^i D. Waters,^^ W.C. Wester lll,^^ 
D. Whiteson^4l A.B. Wicklund,^ E. Wicklund,!^ S. Wilbur,ii F. Wick,^^ H.H. Williams,"^^ J.S. Wilson,36 
P. Wilson,i5 B.L. Winer,^^ R Wittichs,!^ g. Wolbers,!^ H. Wolfe,^^ T. Wright,32 ^ is 2. Wu,^ K. Yamamoto,^^ 
D. Yamato,38 T. Yang,!^ U.K. Yang''," Y.C. Yang,^^ W.-M. Yao,^^ G.P. Yeh,!^ K. Yi",i-"^ J. Yoh,!^ K. Yorita,^'' 
T. Yoshida',38 G.B. Yu,^"^ I. Yu,^^ S.S. Yu,!^ J.C. Yun,!^ A. Zanetti,50 Y. Zeng," C. Zhou,^'^ and S. Zucchelli'=«6 

(CDF CoUaborationt) 

'institute of Physics, Academia Sinica, Taipei, Taiwan 11529, Republic of China 
^Argonne National Laboratory, Argonne, Illinois 60439, USA 
'^University of Athens, 157 71 Athens, Greece 
'^Institut de Fisica d'Altes Energies, ICREA, Uv/iversitat Autonoma de Barcelona, E-08193, Bellaterra (Barcelona), Spain 

^Baylor University, Waco, Texas 76798, USA 
^Istituto Nazionale di Fisica Nucleare Bologna, '^"University of Bologna, 1-40127 Bologna, Italy 
''University of California, Davis, Davis, California 95616, USA 
^University of California, Los Angeles, Los Angeles, California 90024, USA 
^Instituto de Fisica de Cantabria, CSIC- University of Cantabria, 39005 Santander, Spain 
^'^ Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA 
Enrico Fermi Institute, University of Chicago, Chicago, Illinois 60637, USA 
'^Comenius University, 842 48 Bratislava, Slovakia; Institute of Experimental Physics, O4O 01 Kosice, Slovakia 
Joint Institute for Nuclear Research, RU-141980 Dubna, Russia 
'^Duke University, Durham, North Carolina 27708, USA 
'^^^ Fermi National Accelerator Laboratory, Batavia, Illinois 60510, USA 
'^University of Florida, Gainesville, Florida 32611, USA 
'^Laboratori Nazionali di Frascati, Istituto Nazionale di Fisica Nucleare, I-OOO44 Frascati, Italy 
'^University of Geneva, CII-1211 Geneva 4, Switzerland 
'^Glasgow University, Glasgow G12 8QQ, United Kingdom 
'^"Harvard University, Cambridge, Massachusetts 02138, USA 
''''Division of High Energy Physics, Department of Physics, 
University of Helsinki and Helsinki Institute of Physics, FIN-OOOI4, Helsinki, Finland 
University of Illinois, Urbana, Illinois 61801, USA 
'"'"'The Johns Hopkins University, Baltimore, Maryland 21218, USA 
^^Institut fiir Experimentelle Kerriphysik, Karlsruhe Institute of Technology, D-76131 Karlsruhe, Germany 
Center for High Energy Physics: Kyungpook National University, 
Daegu 702-701, Korea; Seoul National University, Seoul 151-742, 
Korea; Sungkyunkwan University, Suwon 440-746, 
Korea; Korea Institute of Science and Technology Information, 
Daejeon 305-806, Korea; Chonnam National University, Gwangju 500-757, 
Korea; Chonbuk National University, Jeonju 561-756, Korea 
Ernest Orlando Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA 
''^'^ University of Liverpool, Liverpool L69 7ZE, United Kingdom 
'''^University College London, London WCIE 6BT, United Kingdom 
Centra de Investigaciones Energeticas Medioambientales y Tecnologicas, E-28040 Madrid, Spain 
^'^Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA 
^'institute of Particle Physics: McCill University, Montreal, Quebec, 
Canada H3A 2T8; Simon Eraser University, Burnaby, British Columbia, 
Canada V5A 1S6; University of Toronto, Toronto, Ontario, 
Canada M5S 1A7; and TRIUMF, Vancouver, British Columbia, Canada V6T 2A3 
University of Michigan, Ann Arbor, Michigan 48109, USA 
Michigan State University, East Lansing, Michigan 48824, USA 
^'^Institution for Theoretical and Experimental Physics, ITEP, Moscow 117259, Russia 
University of New Mexico, Albuquerque, New Mexico 87131, USA 
^'^The Ohio State University, Columbus, Ohio 43210, USA 
Okayama University, Okayama 700-8530, Japan 
Osaka City University, Osaka 588, Japan 
University of Oxford, Oxford 0X1 3RH, United Kingdom 
'^"istituto Nazionale di Fisica Nucleare, Sezione di Padova-Trento, ^ University of Padova, 1-35131 Padova, Italy 



3 



University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA 
'^^Istituto Nazionale di Fisica Nucleare Pisa, University of Pisa, 
'^'^ University of Siena and "Seuola Norrnale Superiore, 1-56127 Pisa, Italy 
'^'^ University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA 
^^Purdue University, West Lafayette, Indiana 41901, USA 
University of Rochester, Rochester, New York 14621, USA 
^^The Rockefeller University, New York, New York 10065, USA 
^"^Istituto Nazionale di Fisica Nucleare, Sezione di Roma 1, 

Sapienza Universita di Roma, 1-00185 Roma, Italy 
^^Rutgers University, Piscataway, New Jersey 08855, USA 
^'^ Texas A&M University, College Station, Texas 11843, USA 

^"istituto Nazionale di Fisica Nucleare THeste/Udine, 
1-34 100 Trieste, ''''University of Udine, 1-33100 Udine, Italy 
''^ University of Tsukuba, Tsukuba, Ibaraki 305, Japan 
•''"^ Tufts University, Medford, Massachusetts 02155, USA 
University of Virginia, Charlottesville, Virginia 22906, USA 

Waseda University, Tokyo 169, Japan 
''^ Wayne State University, Detroit, Michigan 48201, USA 
University of Wisconsin, Madison, Wisconsin 53106, USA 
^^Yale University, New Haven, Connecticut 06520, USA 

Using data from pp collisions at ^/s = 1.96 TeV recorded by the CDF II detector at the Fermilab 
Tevatron, we present improved measurements of the masses and first measurements of natural 
widths of the four bottom baryon resonance states S'^ , S^'^ and , ■ These states are fully 
reconstructed in their decay modes to A^Itv^ where /1° — > ^Jtt" with A'^ pK~n^. The analysis is 
based on a data sample corresponding to an integrated luminosity of 6.0fb~'^ collected by an online 
event selection based on tracks displaced from the pp interaction point. 

PACS numbers: 14.20.Mr, 13.30.Eg, 14.65.Fy 



I. INTRODUCTION 

Baryons with a heavy quark Q as the "nucleus" and 
a light diquark qiq2 as the two orbiting "electrons" can 



* Deceased 

tWith visitors from "Istituto Nazionale di Fisica Nucleare, Sezione 
di Cagliari, 09042 Monserrato (Cagliari), Italy, ''University of CA 
Irvine, Irvine, CA 92697, USA, '^University of CA Santa Barbara, 
Santa Barbara, CA 93106, USA, '^University of CA Santa Cruz, 
Santa Cruz, CA 95064, USA, "^Institute of Piiysics, Academy of 
Sciences of tiie Czech Republic, Czech Republic, ^CERN, CH-1211 
Geneva, Switzerland, ^Cornell University, Ithaca, NY 14853, USA, 
''University of Cyprus, Nicosia CY-1678, Cyprus, 'Office of Sci- 
ence, U.S. Department of Energy, Washington, DC 20585, USA, 
J University College Dublin, Dublin 4, Ireland, '=ETH, 8092 Zurich, 
Switzerland, 'University of Fukui, Fukui City, Fukui Prefecture, 
Japan 910-0017, '"Universidad Iberoamericana, Mexico D.F., Mex- 
ico, "University of Iowa, Iowa City, lA 52242, USA, "Kinki Univer- 
sity, Higashi-Osaica City, Japan 577-8502, ^'Kansas State Univer- 
sity, Manhattan, KS 66506, USA, ^Korea University, Seoul, 136- 
713, Korea, '"UniverswU'Queen Mary, University of London, Lon- 
don, El 4NS, United Kingdom, ^University of Melbourne, Victoria 
3010, Australia, '"Muons, Inc., Batavia, IL 60510, USA, "Nagasaki 
Institute of Applied Science, Nagasaki, Japan, ™ National Research 
Nuclear University, Moscow, Russia, ^Northwestern University, 
Evanston, IL 60208, USA, ^'University of Notre Dame, Notre Dame, 
IN 46556, USA, ^Univcrsidad dc Ovicdo, E-33007 Ovicdo, Spain, 
''''CNRS-IN2P3, Paris, F-75205 France, "''Univcrsidad Autonoma 
de Madrid, Cantoblanco, 28049, Madrid, ''''Texas Tech University, 
Lubbock, TX 79609, USA, <='=Universidad Tecnica Federico Santa 
Maria, llOv Valparaiso, Chile, '''^Yarmouk University, Irbid 211-63, 
Jordan, 



be viewed as the "hehum atoms" of quantum chromody- 
namics (QCD). The heavy quark in the baryon may be 
used as a probe of confinement that allows the study of 
non-perturbative QCD in a different regime from that of 
the light baryons. 

Remarkable achievements in the theory of heavy quark 
hadrons were made when it was realized that a single 
heavy quark Q with mass mq » ^qcd in the heavy 
hadron Hq can be considered as a static color source 
in the hadron's rest frame [1]. Based on this conjecture, 
the light diquark properties of the charm baryon ^4+ (Sc) 
and its bottom partner {^b) can be related by an ap- 
proximate SU(2) symmetry with c ^ b quark exchange. 
Another symmetry emerges because the spin of the heavy 
quark Sq decouples from the gluon field. Models exploit- 
ing these heavy quark symmetries are collectively identi- 
fied as heavy quark effective theories (HQET) [2, 3]. 

As the spin Sqq of a light diquark (plus a gluon field) 
and the spin Sq of a heavy quark are decoupled in HQET, 
heavy baryons can be described by the quantum numbers 
Sq, rriQ, Sqq, m,qq. The total spins of the S'-wave (no or- 
bital excitation) baryon multiplets can be expressed as 
the sum J = Sq + Sqq. Then the singlet yl^ baryon, 
with quark content b[ud] according to HQET, has spin 

of the heavy quark S'f = Its flavor antisymmet- 

ric [ud] diquark has spin Sl^^^ = 0"*" [4]. Under these 
conditions the b quark and the [ud] diquark make the 
lowest-lying singlet ground state = i^. The part- 
ner of the yl^ baryon in the strange quark sector is the 
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yl" baryon. The other two states Ei, and with quark 
content and spin of the flavor symmetric {qq\ diquark 
^{qq} ~ -'^^i constitute two isospin 7=1 triplets with to- 
tal spin = i"*" and — |^ M- These states are the 
lowest-lying S'-wave states that can decay to the singlet 
yl^ via strong processes involving soft pion emission - pro- 
vided sufficient phase space is available. The Eb and El 
particles are classified as bottom baryon resonant states. 
The partners of the states [5] in the strange quark 
sector are E'^*^ baryon resonances, though the = 
E states are light enough to decay only weakly or ra- 
diatively, and only the = |^ states 17(1385) decay 
strongly via the A^-k mode [6] . 

Some recent HQET calculations for bottom baryons 
are available in Ref. [7] . The mass spectra of single heavy 
quark baryons calculated with HQET in combined ex- 
pansions in l/mg and l/A^c, with defined as a num- 
ber of colors, are presented in Ref. [8]. In the potential 
quark model, the mass differences m{SQ) — m{AQ) and 
m{EQ) — m{EQ) are largely due to hyperfine splittings, 
hence the mass differences scale as l/mg. Some recent 
predictions based on potential quark models are found 
in Refs. [9, 10]. There are striking patterns in the masses 
and mass differences of known hadrons. Some of these 
regularities can be understood from known general prop- 
erties of the interactions of quarks, without specifying 
the explicit form of the Hamiltonian. Following this ap- 
proach, the authors of Ref. [11] use semi-empirical mass 
formulae to predict the spectra of c and h baryons. The 
non-perturbativc formalism of QCD sum rules has been 
applied within HQET to calculate the mass spectra of the 
heavy baryons Aq and Eq [12]. Lattice non-relativistic 
QCD calculations for bottom baryons [13] have been 
quite successful, though the uncertainties are typically 
large and exceed the uncertainties of the experimental 
measurements. 

The mass splittings between members of the 1=1 
isospin triplets E^ ' arise from a combination of the in- 
trinsic quark mass difference m{d) > m(?i) and the elec- 
tromagnetic interactions between quarks [10, 14]. Be- 
cause of electromagnetic effects and the d quark be- 
ing heavier than the u quark, the e'^'^ states (with 
composition h\dd\ i.e. all quarks with negative electric 
charge) arc expected to be heavier than the e'^^^ states 
whose composition is h\uu\ [15]. No previous experi- 
mental measurements of isospin mass splitting of bottom 
baryons are available. 

The description of strong decays of baryon resonances 
is a difficult theoretical task [16]. Only a few calcula- 
tions [4, 17, 18] of the E^^ natural widths are available. 
The widths are predicted in the range 4.5 — 13.5MeV/c^ 
for r(Z'6,i"^), and the range 8.5 - 18.0MeV/c2 for 

Until recently, direct observation of h baryons has been 
limited to the yl^ reconstructed in its weak decays to 
J/ipA^ and A'^tt~ [6]. The substantially enlarged experi- 



mental data sets delivered by the Tevatron allow signifi- 
cant advances in the spectroscopy of heavy quark baryon 
states. The resonance Ej^ ' states were discovered by 
CDF [19]. The charged bottom strange S'^ baryon was 
observed and measured [20-22] by both the CDF and 
DO Collaborations. Later, DO reported the first observa- 
tion of the bottom doubly-strange particle Q'^ [23] . Sub- 
sequently the CDF Collaboration confirmed the signal 
and measured the mass of the baryon [22]. Lastly, 
the neutral partner of Sj^, the bottom strange baryon 
H^, was reported for the first time by CDF [24]. Pre- 
cise measurements of the masses and natural widths of 
baryon resonances in the charm sector, specifically the 

Ec , Ec , and yl*"*", were recently reported by the 
CDF Collaboration [25]. 

This study follows the first observation of the Z"^*' 
states using 1.1 fb~^ [19]. We confirm the observation 
of those states using a larger data sample, improve the 

measurement technique, and add new measurements of 

properties of the resonances. In the present analysis, 

the masses of the El '^ and E^ ' states are determined 
independently, with no input from theory assumptions, 
differing from the previous CDF analysis [19]. Using an 
enlarged data sample of 6 fb~^ , we extract the direct mass 
measurements with smaller statistical and systematic un- 
certainties than previously. First measurements of the 
natural widths of the = |^ and = states are 
presented. Based on the new mass measurements, we 
determine the isospin mass splitting for the E^ and E^ 
isospin 7=1 triplets. 

Section II provides a brief description of the CDF II de- 
tector, the online event selection (trigger) important for 
this analysis, and the detector simulation. In Sec. Ill the 
data selection, analysis requirements, and reconstriiction 
of the signal candidates are described. Section IV dis- 
cusses the fit model of the final spectra and summarizes 
the fit results. In Sec. V we estimate the significance of 
signals extracted from the fits. The systematic uncer- 
tainties are discussed in Sec. VI. We present a summary 
of the measurements and conclusions in Sec. VII. 



II. THE CDF II DETECTOR AND SIMULATION 

The component of the CDF II detector [26] most rele- 
vant to this analysis is the charged particle tracking sys- 
tem. The tracking system operates in a uniform axial 
magnetic field of 1.4 T generated by a superconducting 
solenoidal magnet. 

The CDF II detector uses a cylindrical coordinate sys- 
tem with z axis along the nominal proton beam line, ra- 
dius r measured from the beam line and (p defined as an 
azimuthal angle. The transverse plane (r, 4>) is perpen- 
dicular to the z axis. The polar angle, 9, is measured from 
the z axis. The impact parameter of a charged particle 
track do is defined as the distance of closest approach of 
the particle track to the primary vertex in the transverse 
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plane. Transverse momentum, px, is the component of 
the particle's momentum projected onto the transverse 
plane. Pseudorapidity is defined as r] = — ln(tan(^/2)). 

The inner tracking system comprises three silicon 
detectors: layer 00 (LOO), the silicon vertex de- 
tector (SVX II) and the intermediate silicon lay- 
ers (ISL) [27 30]. The innermost part, the LOO detec- 
tor, is a layer of single-sided radiation tolerant silicon 
sensors mounted directly on the beam pipe at a radius 
of 1.35 — 1.6 cm from the proton beam line. It provides 
only an r-0 measurement and enhances the impact pa- 
rameter resolution. Outside this, the five double-sided 
layers of SVX II provide up to 10 track position measure- 
ments. Each of the layers provides an r-(f) measurement, 
while three return a measurement along z, and the other 
two return a measurement along a direction oriented at 
±1.2° to the z axis. The SVX II spans the radii be- 
tween 2.5 cm and 10.6 cm and covers the pseudorapidity 
range |?7| < 2.0. The SVX II detector provides a ver- 
tex resolution of approximately 15 /xm in the transverse 
plane and 70 /xm along the z axis. A fine track impact 
parameter resolution — 35 ^m is achieved, where the 
ado includes an approximate 28 fim contribution from the 
actual transverse size of the beam spot. The outermost 
silicon subdetector, ISL, consists of double-sided layers 
at radii 20 cm to 28 cm, providing two or four hits per 
track depending on the track pseudorapidity within the 
range \r]\ < 2.0 instrumented by the ISL. 

A large open cell cylindrical drift chamber, the central 
outer tracker (COT) [31], completes the CDF detector 
tracking system. The COT consists of 96 sense wire lay- 
ers arranged in 8 superlayers of 12 wires each. Four of 
these superlayers provide axial measurements, and four 
provide stereo views at ±2°. The active volume of the 
COT spans the radial region from 43.4 cm to 132.3 cm. 
The pseudorapidity range \r]\ < 1.0 is covered for tracks 
passing through all layers of the COT, while for the range 
out to 1.0 < |ry| < 2.0, tracks pass through less than 
the full 96 layers. The trajectory of COT tracks is ex- 
trapolated into the SVX II detector, and the tracks are 
refitted with additional silicon hits consistent with the 
track extrapolation. The two additional layers of the ISL 
help to link tracks in the COT to hits in the SVX II. 
The combined track transverse momentum resolution is 
(t{pt)/pt ^ 0.07% pt [GeY/c]-K 

The analysis presented here is based on events recorded 
with a three-tiered trigger system configured to collect 
large data samples of heavy hadrons decaying through 
multi-body hadronic channels. We refer to this as the 
displaced two-track trigger. We use two configurations 
of this trigger, the "low-px" and the "medium-pT" selec- 
tions. At level 1, the trigger uses information from the 
hardware extremely fast tracker [32] . The "low-px" con- 
figuration of the displaced two-track trigger requires two 
tracks in the COT with pt > 2.0 GeV/c for each track, 
and with an opening angle of \A(j)\ < 90° between the 
tracks in the transverse plane. Additionally the track 
pair scalar sum must satisfy pxi +PT2 > 4-0 GeV/c. The 



corresponding criteria imposed in the "medium-px" con- 
figuration are px > 2.0 GeV/c for each track, opening 
angle |A</)| < 135°, and pxi + PT2 > 5.5GeV/c. The 
level 2 silicon vertex trigger (SVT) [33, 34] associates the 
track pair from the extremely fast tracker with hits in 
the SVX II detector and recognizes both tracks using a 
large look-up table of hit patterns. The SVT repeats 
the level 1 pt criteria and limits the opening angle to 
2° < \A(p\ < 90°. Only in the case of the medimn-j>r 
configuration are the charges of the tracks required to be 
of opposite sign. Crucially, the SVT imposes a require- 
ment on the transverse impact parameter of each track 
to be 0.12 < do < 1mm, given the excellent resolution 
provided by SVX II. Finally, the distance in the trans- 
verse plane between the beam axis and the intersection 
point of the two tracks projected onto their total trans- 
verse momentum is required to be L^y > 200 /xm. The 
level 3 software trigger uses a full reconstruction of the 
event with all detector information and confirms the cri- 
teria applied at level 2. The trigger criteria applied to 
the do of each track in the pair and to L^y preferentially 
select decays of long-lived heavy hadrons over prompt 
background, ensuring that the data sample is enriched 
with b hadrons. 

The mass resolution on the ' resonances is pre- 
dicted with a Monte Carlo simulation that generates 
b quarks according to a next-to-leading order calcu- 
lation [35] and produces events containing final state 
hadrons by simulating b quark fragmentation [36]. Mass 
values of 5807.8 MeV/c^ for Eb and 5829.0 MeV/c^ for 
[19] are used in the Monte Carlo generator. Final 
state decay processes are simulated with the evtgen [37] 
program, and all simulated b hadrons are produced with- 
out polarization. The generated events arc input to the 
detector and trigger simulation based on GEANt3 [38] 
and processed through the same reconstruction and anal- 
ysis algorithms as are used on the data. 



III. DATA SAMPLE AND EVENT SELECTION 

This analysis is based on data equivalent to 6.0 fb^^ of 
pp collisions collected with the displaced two-track trigger 
between March 2002 and February 2010. We study S^,*^ 

resonances in the exclusive strong decay mode Z'^*-'^ 
A'f^nf, where the low momentum pion is produced 
near kinematic threshold [39]. The A'^ decays to A'^tt^ 
with a prompt pion tt^ produced in the weak decay. This 
is followed by the weak decay yl+ — > pK~ Tr~^. 

To reconstruct the parent baryons, the tracks of 
charged particles are combined in a kinematic fit to form 
candidates. No particle identification is used in this anal- 
ysis. The following two complementary quantities de- 
fined in the plane transverse to the beam line and relat- 
ing the decay path of baryons to their points of origin are 
used: the proper decay time of the baryon candidate h 
expressed in length units ct{h), and the impact parame- 
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ter do{h). Specifically, the decay length is defined as 



ct{h) = Lxy{h) 



M(h) c 
PT{h) 



(1) 



where L^yih) is expressed in length units and defined 
as the projection onto PT{h) of the vector connecting 
the primary vertex to the heavy baryon decay vertex in 
the transverse plane. The transverse impact parameter 
do{h) of the candidate is defined analogous to the one 
of a charged particle track. An event-specific primary 
interaction vertex is used in the calculation of the ct{h) 
and d(i{h) quantities. The measurement uncertainties Get 
and (Tda originate from the track parameter uncertainties 
and the uncertainty on the primary vertex. 



TABLE I: Analysis requirements for ■^t'^b reconstruc- 

tion. The quantity ct{Af <— A^) is defined analogously to 
Eq. (1) as the proper time where Lxy{A'^) is calculated 
with respect to the A° vertex. 



Quantity 



Requirement 



> 200 nm 

> 12.0 

< 80 /xm 

> -150 /um 

< 250 fim 

PT(vr-) > 1.5 GeV/c 

PtK') >4.0GeV/c 
Prob(x3i3) of A^ vertex fit > 0.01% 



ctiA'i)/act 
do{At) 
ct{A+ ^ Al) 
ct{At ^ A°b) 



A. Reconstruction of the A^ candidates 

The analysis begins with reconstruction of the yl+ — >■ 
pK~Tr~^ decay by fitting three tracks to a common vertex. 
The invariant mass of the /l^ candidate is required to 
be within ±18MeV/c^ of the world-average yl+ mass [6]. 
The momentum vector of the yl+ candidate is then ex- 
trapolated to intersect with a fourth pion track, the 
TT^-candidate, to form the A'^ — > A^tt^ candidate ver- 
tex. The A'-^ vertex is subjected to a three-dimensional 
kinematic fit with the yl+ candidate mass constrained 
to its world average value [6]. The probability of the 
constrained yl^ vertex fit must exceed 0.01%. Standard 
quality requirements are applied to each track, and only 
tracks with pt > 400 MeV/c are used. All tracks are re- 
fitted using pion, kaon and proton mass hypotheses to 
properly correct for the differences in multiple scattering 
and ionization energy loss. At least two tracks among 
the p, K~, TT"*", and tt^ candidates are required to fulfill 
the level 2 (SVT) trigger requirements. 

To suppress prompt backgrounds from the primary in- 
teraction, the decay vertex of the A° is required to be 
distinct from the primary vertex. To achieve this, cuts 
on cf(yl{J) and its significance ci(/l°)/(Tci arc applied. We 
require the A'^ vertex to be close to the A*^ vertex by 
applying cuts on ct{A^) where the corresponding quan- 
tity Lxy(A'^) is calculated with respect to the /l^ ver- 
tex. The requirement ct{A'^) > —150 /xm reduces con- 
tributions from yl+ baryons directly produced in pp in- 
teraction and from random combination of tracks faking 
yl+ candidates which may have negative ci(/l+) values. 
The other restriction, ct(yl+) < 250 /im, aims at reduc- 
ing contributions from B° — >• D+tt" decays, followed by 
— iC~7r+7r+ decays. The requirements take into ac- 
count ct resolution effects and exploit the much shorter 
A^ lifetime compared to the £>+ [19, 40]. To reduce 
combinatorial background and contributions from par- 
tially reconstructed decays, we ask A'^ candidates to point 
to the primary vertex by requiring the impact parame- 
ter do{A^) not to exceed 80 /um. The choice of analy- 
sis requirements to identify yl° — > A'^tt^ candidates is 
made using an optimization based on the experimental 



data only. The figure of merit S/\/S + B is used dur- 
ing the optimization, where S is the A^ signal and B is 
the background under the signal, respectively. At every 
step of the optimization procedure, both quantities are 
obtained from fits of the A'^n^ invariant mass spectrum 
and are determined from the corresponding numbers of 
candidates fit within ±3a of the /l" signal peak. Table I 
summarizes the resulting /Ijj analysis requirements. 

Figure 1 shows a prominent A'^ signal in the A'^tt'^ 
invariant mass distribution, reconstructed using the op- 
timized criteria. A binned maximum-likelihood fit finds 
a signal of approximately 16 300 candidates at the ex- 
pected /l^ mass, with a signal to background ratio around 
1.8. The fit model describing the invariant mass distri- 
bution comprises the Gaussian A^ A'^tt^ signal on 
top of a background shaped by several contributions. 
Random four-track combinations dominating the right 
sideband are modeled with an exponentially decreasing 
function. Coherent sources populate the left sideband 
and leak under the signal. These include reconstructed 
B mesons that pass the yl" — )• A^n^ selection crite- 
ria, partially reconstructed A^ decays, and fully recon- 
structed A^ decays other than A'^tt^ [e.g. yl" A'^K^). 
Shapes representing the physical background sources are 
derived from Monte Carlo simulations. Their normaliza- 
tions are constrained to branching ratios that are either 
measured (for B meson decays, reconstructed within the 
same A^tt,^ sample) or theoretically predicted (for A^ 
decays) [19, 40]. 



B. Reconstruction of Ej^*' candidates 

To reconstruct the l^^*-*^ A'^irf candidates, each 
A^TT^ candidate with invariant mass within the A'^ sig- 
nal region, 5.561 — 5.677 GeV/c^, is combined with one 
of the tracks remaining in the event with transverse mo- 
mentum down to 200 McV/c. The A'^ mass range covers 
±3 standard deviations as determined by a fit to the 
signal peak of Fig. 1. To increase the efficiency for recon- 
structing Sj^ decays near the kinematic threshold, the 
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FIG. 1: Invariant mass distribution of ylf, — >■ A^TTf^ candidates 
with the projection of a mass fit overlaid. 



TABLE II: X"^* candidate selection requirements. 

Quantity Requirement 

m{AtTT^) G (5.561, 5.677) GeV/c^ 

doi'T't) < 0.1 cm 

PT(7r±) > 200 MeV/c 

dQ{nt)/ado < 3.0 

PT{Tvt) <Pt(V) 

PTjSl*^^) > 4.0 GeV/e 



quality criteria applied to soft pion tracks are loosened in 
comparison with tracks used for the A*^ candidates. The 
basic COT and SVX II hit requirements are imposed on 
TT^ tracks, and only tracks with a valid track fit and error 
matrix are accepted. 

Random combinations of signal candidates with 
tracks constitute the dominant background to the 
Z"^*-*^ ^ '^b''^t signal. The remaining backgrounds are 
random combinations of soft tracks with B mesons recon- 
structed as A*^ baryons, and combinatorial background 
events [19]. To reduce the background level, a kine- 
matic fit is applied to the resulting combinations of A'^ 
candidates and soft pion tracks to constrain them 
to originate from a common point. Furthermore, since 
the bottom baryon resonance originates and decays at 
the primary vertex, the soft pion track is required to 
point back to the primary vertex by requiring an impact 
parameter significance, do{Trf)/adQ, smaller than three. 
The transverse momentum of the soft pion is required to 
be smaller than the tt^ transverse momentum. As we 
already require (tt^ ) > 1.5GeV/c (Table I) the condi- 
tion imposed on the soft pion is fully efficient. The 

Z'^ candidate selection requirements are summarized 
in Table II. 



IV. DETERMINATION OF RESONANCE 
PROPERTIES 

The analysis of the mass distributions is per- 

formed using the Q value 

Q = miA",7rt) - m(^g) - , (2) 

where TOjt is the known charged pion mass [6] and rn(yl^) 
is the reconstructed A'^n^ mass. The mass resolution of 
the A^ signal and most of the systematic uncertainties 
cancel in the mass difference spectrum. The and 
signals are reconstructed as two narrow structures 
in the Q-value spectrum. The properties, yields, and 
significance of the resonance candidates are obtained by 
performing unbinned maximum-likelihood fits on the Q- 
value spectra. 

The shapes of the Z'^ resonances are each modeled 
with a non-relativistic Breit-Wigner function. Since the 
soft pion in Z"^ strong decay modes is emitted in a P- 
wave, the width of the Breit-Wigner function is modified 
as follows [41]: 

r(Q;Qo,ro) = ro , (3) 

where Qo is the Q value at the resonance pole; p* ^ and "p*^ 

are the momenta of the soft pion in the Z^ rest frame, 
off and on the resonance pole respectively; and Fq is the 
corrected width. The soft pion momenta are calculated 
based on two-body decay kinematics [6]. Both Qo and 
Fq are floating flt parameters. 

The Breit-Wigner function is convoluted with the de- 
tector resolution, which is described by a narrow core 
Gaussian plus a broad Gaussian. Their widths f7„ and 
and relative weights g„ and (1 — g^) are calculated 
from the CDF full Monte Carlo simulation. Numerical 
convolution is necessary because the modified width de- 
pends on the mass. The effects of imperfect modeling in 
the simulation are discussed with the systematic uncer- 
tainties in Sec. VI. 

We use a kinematically motivated model for the back- 
ground, described by a second order polynomial modu- 
lated with a threshold square root-like term. 



r\Q;C,bi,b2), 

where C, bi, and 62 are the second order polynomial 
coefficients and mx is a threshold fixed to 0.140 GeV/c^, 
the mass of the pion. 

The full model for the Q-value spectra of all isospin 
partner states Z^ and Z^ describes two narrow 
structures on top of a smooth background with a thresh- 
old. The negative logarithm of the extended likelihood 
function (NLL) is minimized over the unbinned set of Q 
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values observed for N candidates in data: 

N 

- In (£) = - ^ ln(Ari Si + N2S2 + Nb BQ) 

(5) 

+ (TVi + N2 + m) 

-Nln{Ni+N2+Nh). 

Independent likelihood functions are used for 17^*^'*' and 

Ej^ ' candidates. The Q- value spectrum is fit over the 
range 0.003 - 0.210 GeV/c^. The effect of this choice is 
discussed in Sec. VI. The probability density functions 
(PDF) in Eq. (5) are defined as follows: 

{i) Si = S{Q;Ql,rl,ai,gi,ai,) is the normal- 
ized convolution of a Breit-Wigncr and a double 
Gaussian responsible for the i^b ) (* ^ ^) 
S^'^ {i = 2) signals. Here Qq is the floating 

pole mass and Fq is the floating natural width. The 
detector's Gaussian resolution parameters cr^, crjy 
and 5^ are set from the Monte Carlo data. A dom- 
inant with Qn ~ 70% relative weight narrow core 
an of about 1.2MeV/c^ is set for the (Ej^) and 
about 1.4MeV/c2 for S*+ {S*-). A broad compo- 
nent a-u, of about 2.9MeV/c^ is set for the {^b) 
and about S.SMeV/c^ for r*+ (^T)- 

(ii) Ni is the floating yield of the {U^} {i = 1) or 

s;+ (u;-) (» = 2). 

(iii) Bg = Bg{Q;mT,CMM) is the PDF correspond- 
ing to the background form in Eq. (4). 

(iv) Nb is the floating yield of the background contri- 
bution. The sum of fitted yields, Ni -\- N2 + Nb, 
is the Poisson mean value of the total number of 
candidates N for the particular species , or 
S^, corresponding to isospin triplets Sb and 

The total number of floating parameters in the fit per 
each pair of isospin partners is nine. 

Extensive tests on several thousand statistical trials 
show that the likelihood fit yields unbiased estimates 

with proper uncertainties. 

The experimental sj^*^ and U^*^^ Q-value distribu- 
tions, each fitted with the unbinned likelihoods described 
above, are shown in Fig. 2. The projection of the cor- 
responding likelihood fit is superimposed on each graph. 
The Q- value distributions show clear signals of , 
and S'^ , -^6^) respectively. The pull distributions are 
shown in the bottom plots of both figures and are cal- 
culated as the residuals of the histogram with respect to 
the corresponding likelihood fit projection normalized by 
the data uncertainty. Both pull distributions are evenly 
distributed around zero with fluctuations of ±2(t, approx- 
imately. The fit results are given in Table III. 



TABLE III; Summary of the results of the fits to the 
Q = M(ylj7r*) - M(Al) - m-n spectra. The statistical un- 
certainties are returned by the unbinned maximum-likelihood 
fits. 



State Qo value, Natural width. Yield 
MeV/c^ To, MeV/c^ 





56.2±;;;« 


4 9+a-i 
^■^-2.1 


340+^,^ 




75.8 ±0.6 


1 -0-1.8 


5401^° 




ro 1 +0.9 
'-'^•-'--0.8 


9.71-1 


4701-° 


5-,*+ 


72.8 ± 0.7 


11.511^ 


800llJ° 



V. SIGNAL SIGNIFICANCE 

The significance of the signals is determined using a 
log-likelihood ratio statistic [42, 43], -2 \a{Co/Li) . We 
define hypothesis Hi corresponding to the presence of 
S'^ , El~ or S'^ , El'^ signals on top of the background. 
The 'Hi hypothesis is described by the likelihood C\; see 
Eq. (5). The various null hypotheses, each identified with 
"Ho and nested to Hi correspond to a few different less 
complex scenarios described by the likelihood jCq. The 
likelihood ratio is used as a variable to derive p val- 
ues for observing a deviation as large as is in our data 
or larger, assuming Ho is true. The number of degrees 
of freedom of the equals the difference AA^^of in the 
number degrees of freedom between the Hi and Ho hy- 
potheses in each case. We consider the following types 
of Ho to estimate the significance of the two-peak signal 
structure and of individual peaks of the observed Sj^ ' 
and U^*'"*" states: 

(i) A single enhancement is observed anywhere in the 
fit range. The corresponding likelihood £0 includes 
only a single peak PDF on top of the background 
form in Eq. (4), the same as for the Ci. The 
difference in the number of degrees of freedom is 
AA^dof = 3. The width Fq floats in the fit over the 
wide range 1 — 70MeV/c^. The position of the en- 
hancement Qo is allowed to be anywhere within the 
default fit range. We test the case in which the ob- 
served two narrow structures could be an artifact of 
a wide bump where a few bins fiuctuated down to 
the background level. 

(ii) The signal is observed but the Sb is interpreted 
as background. We impose a loose requirement on 
the existence of the second peak, E'^ fixing only 
the width of E^ to the expected theoretical value 
of 12MeV/c2 [17]. We let the fitter find the 
position within the default fit range. The number 
of free parameters is changed by 4. 

(iii) The signal Eb is observed but the E^ is interpreted 
as background. This null hypothesis is similar to 
the previous one. The width of the Eb is fixed to 
7MeV/c2 [17]. 
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FIG. 2: The left (right) plot shows the Q-value spectrum for S^*^ (X'^*'^) candidates with the projection of the corresponding 
unbinned likelihood fit superimposed. The Q value is defined in Eq. (2). The pull distribution of each fit is shown in the bottom 
of the corresponding plot. 



(iv) Neither the 17^ nor the SI is observed, and the "Hq 
hypothesis is the default background model used in 
£.1. We consider the case in which the smooth back- 
ground fluctuates to two narrow structures corre- 
sponding to the T-Li hypothesis. The difference in 
the number of degrees of freedom is 6. 

In addition to all the cases considered above, we intro- 
duce an additional case in which the Hi hypothesis cor- 
responds to any single wide enhancement considered in 
(i) while the "Ho hypothesis is the default background 
considered in (iv). This special test determines the sig- 
nificance of the single enhancement with respect to pure 
background. 

Table IV summarizes the results of these tests. The 
null hypothesis most likely to resemble our signal is a 
broad single enhancement fluctuating to the two narrow 
structures. The results of this study establish conclu- 
sively the and ' signals with significance of 

6a or higher. 



TABLE IV: Statistical significances of the observed signals 
against various null hypotheses. Na is the calculated number 
of Gaussian standard deviations based on Prob(x^). 



Wo 


States 




Hi 


Any single wide 




6.7 


Two narrow 


enhancement 




6.1 


structures 


No structures 


^b 


10.7 Any single wide 




^b 


13.2 


enhancement 


No X'i,, with X'j, 


y^(*)- 
^b 


7.6 


Two narrow 


ro2 = 12 MeV/c^ 


^b 


7.9 


structures 


No El, with Eb, 


^b 


10.0 


Two narrow 


Toi = 7MeV/c2 


^b 


12.5 


structures 


No structures 


^b 


12.4 


Two narrow 




^b 


14.3 


structures 



VI. SYSTEMATIC UNCERTAINTIES 

The systematic uncertainties considered in our analysis 
are the following: 

(i) The uncertainty due to the CDF tracker momentum 
scale. 

(ii) The uncertainty due to the resolution model (see 
Sec. IV) described by the sum of two Gaussians. 



This source is expected to dominate the systematic 
uncertainties on width measurements. 

(iii) The choice of background model. 

(iv) An uncertainty due to the choice of Q-value fit 
range. 

To calibrate the tracker momentum scale, the energy 
loss in the material of CDF tracking detectors and the 
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strength of the magnetic field must be determined. Both 
effects are calibrated and analyzed in detail using high 
statistics samples of J/ip , ^'(25'), T{1S), reconstructed 
in their decay modes as well as — >• K~'it^, 

tp{2S) J/ip{^ ^i+n~)-K+-K~ [44, 45]. The correspond- 
ing corrections are taken into account by tracking algo- 
rithms. Any systematic uncertainties on these correc- 
tions are largely negligible in the Z"^ Q-value measure- 
ments. The uncertainties on the measured mass differ- 
ences due to the momentum scale are estimated from 
the deviations between Qq values observed in similar 
decays reconstructed in CDF data and the known Qo 
values [6]. The reference modes are S'^^ — >• vl^7r+, 
i:° ^ A+-K-, A*+ yl+TT+TTj, and D*+ D^-Kt . The 
linear extrapolation of the measured offsets as a function 

of Qo towards the Z"^ kinematic regime is taken as the 
mass-scale uncertainty. The determined systematic un- 
certainty on the momentum scale covers also any residual 
charge-dependence of the scale. For the mass difference 
Qo, the systematic uncertainty due to a possible imper- 
fect alignment of the detector is negligible [44] . 

Following the method used in Ref. [46], the D*+ 
D^{-^ ii'~7r+)7r+ signal peak in the mass difference dis- 
tribution m{D*'^) — m{D^) has been reconstructed in sev- 
eral bins of soft pion transverse momentum pt{iTs) start- 
ing with 200MeV/c as in the data. Each signal distri- 
bution is subjected to an unbinned maximum-likelihood 
fit with the sum of a Breit-Wigner function convoluted 
with a double Gaussian function to describe the detec- 
tor resolution. The background under the D*^ signals 
is described by an empirical function [47, 48]. For each 
of the P't{'Ks) bins, the fit determines the D*^ width, 
which never exceeds 0.2MeV/c^. Because the D*^ nat- 
ural width is much smaller than the tracking resolution, 
the value of 0.2MeV/c^ is assigned as a systematic un- 

certainty on the measured ' natural width due to the 
momentum scale of the CDF tracker. 

Unless otherwise specified, the systematic uncertain- 
ties discussed below are evaluated for the measurable 
quantities Qo and Fq by generation of statistical trials. 
In each trial, the sample is generated according to the 
PDF (see Table III) with the nuisance parameters mod- 
ified by the uncertainty with respect to the default set 
of parameters. Then the sample is subjected to the un- 
binned maximum-likelihood fit twice, with the default 
PDF and with the PDF of the modified nuisance param- 
eter set. The fit results are compared on a trial-by-trial 
basis, and their difference is computed. The systematic 
uncertainty is found from the mean of a Gaussian fit of 
the distribution of the computed differences. 

The statistical uncertainties on the resolution model 
parameters due to the finite size of the Monte Carlo 
datasets introduce a systematic uncertainty. Variations 
of the double Gaussian widths C7„ and Uw and the weight 
Qn within their statistical uncertainties returned from the 
fits of Monte Carlo spectra are propagated into the mea- 
surable quantities using the statistical trials. 



The CDF tracking simulation does not reproduce with 
perfect accuracy the tracking resolutions, especially for 
soft tracks at the kinematic threshold of decays. 
To estimate this contribution, we use the D*^ me- 
son decay as the reference mode reconstructed down to 
PT{-Kf) ~ 200MeV/c in the observed and simulated sam- 
ples. We compare the mass resolution of the reference 
signal found in data with the one predicted by Monte 
Carlo simulation. The comparison is made independently 
for D*+ -D°7r+ and D*^ D^t^J states, as a function 
of soft pion pt using early data (Period 1) and late data 
(Period 2). Figure 3 shows the comparisons of the narrow 
core resolution between the data and Monte Carlo both 
for D*+ (left plot) and £>*" (right plot). The resolution 
is stable as a function of data-taking time. 

The CDF Monte Carlo simulation typically under- 
estimates the D*^ resolutions in the experimental 
data: (T„(data) < 1.25 (t„ (Monte Carlo). Similar rela- 
tions are found for the broad component of the resolu- 
tion: ^^(data) < 1.40 (T,„( Monte Carlo). These factors 
are used as the sources of the systematic uncertainties. 
The resolution extracted for the D*~ is systematically 
smaller than for the £)*+ by at most 20% for o",i and by 
at most 40% for a^. The Monte Carlo predictions for a„ 
and CTu, are decreased by these latter factors to estimate 
the other bounds of the systematic miccrtainties. In both 
cases the conservative approach is taken. 

To find the systematic uncertainty associated with the 
choice of background shape, we change our background 
PDF to the one used for the D*^ mass difference spec- 
tra [47, 48] and compare with the default background 
PDF. 

The uncertainty associated with the fit range is es- 
timated by varying the default low edge down to 
0.0015 GcV/c2 and up to 0.006 GeV/c^. The fit results 
are slightly sensitive to the choice of the low edge and 
any observed biases are assigned as another systematic 
uncertainty. 

The final systematic uncertainties are listed in Table V. 



VII. RESULTS AND CONCLUSIONS 

The analysis results are arranged in Table VI. Prom 
the measured Ef^ Q values we extract the absolute 
masses using the known value of the tt^ mass [6] and 
the CDF mass measurement, m(ylg) = 5619.7 ± 
1.2 (stat) ±1.2 (syst) McV/c^ , as obtained in an indepen- 
dent sample [44]. The yl^ statistical and systematic un- 
certainties contribute to the systematic uncertainty on 
the Ej^ ' absolute masses. 

Using the measured Q values, we extract the isospin 
mass splittings for the isotriplets of the = and 

= 1^ states. The statistical uncertainties on the 

Q-measurements of the corresponding charge states are 
added in quadrature. We assume that the correlated sys- 
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FIG. 3: The left (right) plot shows the ratio of the widths of the narrow component of the D*"*" (D*~) mass resolution for data 
and simulation (circles) and for different subsamples of data (triangles) as a function of the transverse momentum of the soft 
pion. The last bin on every plot corresponds to a statistics integrated above l.OGeV/c. 



TABLE V: Summary of the systematic uncertainties listed in the following order: mass scale, resolution, choice of background 
model, and fit range. The total systematic uncertainty is obtained by adding all the associated uncertainties in quadrature. 
The last column shows the percentage of the total systematic uncertainty relative to its central value. 



Measurable quantity 


Scale 


Resolution Background Fit range 


Total 


Percentage 


Q(r-)[MeV/c2] 


-0.38 


+0.06 
-0.07 


+0.04 
-0.04 


+0.02 
-0.03 


+0.07 
-0.39 


+0.1 
-0.7 


TiU-)lMeY/c'] 


+0.20 
-0.20 


+0.85 
-0.87 


+0.50 
-0.50 


+0.50 
-0.51 


+1.13 
-1.14 


+23 
-23 


Q(r*-)[MeV/c2] 


-0.56 


+0.06 
-0.08 


+0.06 
-0.06 


+0.02 
-0.09 


+0.09 
-0.58 


+0.1 
-0.8 


r(i;r)[MeV/c2] 


+0.20 
-0.20 


+0.65 
-0.96 


+0.30 
-0.30 


+0.50 
-0.90 


+0.89 
-1.36 


+ 12 

-18 


Q(r+)[MeV/c2] 


-0.35 


+0.07 
-0.12 


+0.05 
-0.05 


+0.02 
-0.03 


+0.09 
-0.38 


+0.2 
-0.7 


r(r+)[MeV/c2] 


+0.20 
-0.20 


+0.94 
-0.90 


+0.40 
-0.40 


+0.50 
-0.51 


+1.16 
-1.12 


+ 12 
-12 


Q{S;+)[MeV/c^] 


-0.52 


+0.06 
-0.13 


+0.10 
-0.10 


+0.02 
-0.09 


+0.12 
-0.55 


+0.2 
-0.8 


r(r*+)[MeV/c2] 


+0.20 
-0.20 


+0.64 
-1.01 


+0.50 
-0.50 


+0.50 
-0.90 


+0.97 
-1.46 


+8.5 
-13 



tematic uncertainties due to mass scale, fit bias due to 
choice of fit range, and imperfect Monte Carlo description 
of the resolution are completely canceled in the isospin 
mass splittings. The uncertainties due to background 
choice are added in quadrature. 

In conclusion, we have measured the masses and widths 
of the baryons using a sample of approximately 

16 300 A'^ candidates reconstructed in their A^t^^ 
mode corresponding to 6fb^^ of CDF data. 

The first observation [19] of the s'^^^ bottom baryons 
has been confirmed with every individual signal recon- 
structed with a significance well in excess of six Gaussian 
standard deviations. 

The statistical precision on the direct mass differences 



is improved by a factor of two over the previous mea- 
surement [19]. The measurements are in good agreement 
with the previous results and supersede them. 

The isospin mass splittings within the 7 = 1 triplets 

of the Sb and states have been extracted for the first 

f *) — 

time. The states have higher masses than their 

Sf^ partners, following a pattern common to most of 
the known isospin multiplets [15]. This measurement fa- 
vors the phenomenological explanation of this ordering 
as due to the higher masses of the d quark with respect 
to the u quark and the larger electromagnetic contribu- 
tion due to electrostatic Coulomb forces between quarks 
in SI states than in Sf^ ones. The difference in the 
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TABLE VI: Summary of the final results. The first uncertainty is statistical and the second is systematic. 



State Q value, MeV/c^ Absolute mass m, MeV/c^ Natural width T, MeV/c^ 



^6 


n +0.6 +0.1 

00-'' -0.5-0.4 


5815.5 Zo.5 ± 1-7 


d n -l-3. 1 1 -1 -1 

4.9 T2.1 =t 1.1 




i-TC o 1 r\ -|-0.1 

75.8 ± 0.6 _Q g 


5835.1 ± O.o _j g 


7 c +2.2+0.9 
' -0 -1.8 -1.4 




C-r, 1 +0.9+0.1 

OZ.i -0.8 -0.4 


c:qi 1 q +0.9 i -t n 
0811. _o 8 ± ' 


Q 7+3.8+1.2 
' -2.8 -1.1 




72.8 ±0.7lo;6 


5832.1 ±0.7li;8 


r+2. 7+1.0 
ll.O -2.2 -1.5 






Isospin mass splitting, 


MeV/c^ 


m{E+) - m{E^) 




-4.2t{;J±0.1 




m{s;+)-mis;-) 




-3.0lJ:°±0.1 





measured isospin mass splittings between the and 

isotriplets supports the theoretical estimate of Ref. [10]. 
The natural widths of the and states have been 
measured for the first time. The measurements are in 
agreement with theoretical expectations. 
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